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Artificially confined, small quantum systems show a 
high potential for employing quantum physics in tech- 
nology. Ultra-cold atom gases have opened an excit- 
ing laboratory in which to explore many-particle sys- 
tems that are not accessible in conventional atomic or 
solid state physics 1 . It appears promising that optical 
trapping of cold bosonic or fermionic atoms will make 
construction of devices with unprecedented precision 
possible in the future, thereby allowing experimenters 
to make their samples much more "clean" , and hence 
more coherent. Trapped atomic quantum gases may 
thus provide an interesting alternative to the quan- 
tum dot nanostructures 2 produced today. Optical 
lattices created by standing laser waves loaded with 
ultra-cold atoms are an example 1 ' 3 of this. They pro- 
vide a unique experimental setup to study artificial 
crystal structures with tunable physical parameters. 
Here we demonstrate that a two-dimensional opti- 
cal lattice loaded with repulsive, contact-interacting 
fermions shows a rich and systematic magnetic phase 
diagram. Trapping a few (TV < 12) fermions in each of 
the single-site minima of the optical lattice, we find 
that the shell structure in these quantum wells deter- 
mines the magnetism. In a shallow lattice, the tunnel- 
ing between the lattice sites is strong, and the lattice 
is non-magnetic. For deeper lattices, however, the 
shell-filling of the single wells with fermionic atoms 
determines the magnetism. As a consequence of 
Hund's first rule, the interaction energy is lowered by 
maximizing the number of atoms of the same species. 
This leads to a systematic sequence of non-magnetic, 
ferromagnetic and antiferromagnetic phases. 
An optical lattice resembles an "egg box" -like arrange- 
ment of single quantum wells, each confining a small 
number of atoms. Different lattice geometries can be 
realized (see D. Jaksch and P. Zoller— for a review). The 
tunneling and the localization of atoms in the lattice are 
controlled by the lattice depth which can be tuned by 
changing the laser intensity. This allows for a smooth 
transition from a tightly bound lattice to a system of 
nearly free atoms. The confined atoms have many inter- 
nal (hyperfine) states which can be manipulated by laser 
light. In a weakly interacting, dilute atom gas, s-wave 
scattering dominates. The strength of the short-ranged 
contact interaction between the atoms can be tuned in 
the vicinity of a Feshbach resonance 4,5 ' 6 ' 7 ' 8,9 . 



For bosons in an optical lattice, it was possible to 
realize the Mott insulator-superfluid quantum phase 
transitionifiiiiiSiii More recently, a fermionic many- 
particle quantum system on a three-dimensional lattice 
was experimentally studied by Kohl et a/J^l 5 -, who inves- 
tigated the transition from a band insulator to a normal 
state. 

The atom dynamics in an optical lattice is often described 
by the Hubbard or the Bose- Hubbard models 1 ^ with the 
hopping parameter t and the on-site interaction parame- 
ter U . These models rely on the single-band approxima- 
tion which allows at most two atoms in a lattice site in the 
fermionic case. From an experimental viewpoint, it has 
been argued that these discrete models are ideal for de- 
scribing contact-interacting atoms in an optical lattice 10 . 
In condensed matter physics, the Hubbard model is con- 
stantly used to describe magnetic correlations in solids. 
In a deep lattice (with small t/U) for a half-filled band, or 
one fermion per site, the Hubbard model predicts antifer- 
romagnetism. In fact, it can be shown that in this case 
the Hubbard Hamiltonian coincides with the antiferro- 
magnetic Heisenberg model of localized atoms 1 *. When 
the band is filled, each lattice site carries two fermions 
with opposite spins, and magnetism can not be observed. 
More recently, a SU(n) Hubbard model was applied to de- 
scribe fermionic atoms with n different "flavors" or spin 
states 1 ^. 

Here, we go beyond the lowest-band approximation, ap- 
plying the mean- field approach for a two- component sys- 
tem: We consider fermionic atoms with two hyperfine, 
or spin, species confined into a two-dimensional (square) 
lattice, where it is possible to trap a few (N < 12) atoms 
in each single well (at each single site). We find that 
depending on the statistics and the spin of the trapped 
atoms, the optical lattices can have an intriguing and rich 
magnetic structure. The magnetism follows closely the 
shell structure in the individual lattice wells, emerging as 
there is an excess of one spin species in a single lattice 
site. The situation is in fact analogous to the behavior of 
quantum dot lattices in semiconductors, where the mag- 
netism is governed by the shell structure of the single-dot 
comp onent si 9 ^ ^!. 

An optical lattice can be created by counter-propagating 
laser beams. The resulting potential is a sinusoidal stand- 
ing wave V opt (r) — Vb(cos 2 (fccc) + cos 2 (fey) ), where the 
amplitude Vq is tuned by laser intensity and the wave 
number k = 2tt/X is set by the laser wave lengthii. A 
natural unit for the energy is the recoil energy Er = 
H 2 k 2 /(2m), and length is measured by the inverse of 
the wave number. (In experiments, an underlying slowly 
varying harmonic confinement is added giving rise to an 
inhomogeneous filling of the lattice and coexistence of 
insulating and conducting domains. This external con- 
finement can be made small and is neglected here.) 
We consider a square lattice with two lattice sites in the 
unit cell, as sketched in Figure 1(A) and (B). One site 
resides at the center of the unit cell, and the other one 
crosses the corner periodically, as indicated in the con- 
tour plot of the optical potential. Note that this is the 
simplest choice of the unit cell allowing for antiferromag- 
netic alignment of the single-site spins. The inter-site 
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FIG. 1: (A) Schematic sketch of the optical lattice potential, 
V pt, confining a small number of atoms in each of the po- 
tential minima. (Here, the arrows indicate a ferromagnetic 
ground state as an example). (B) shows the contours of V op t 
and indicates the unit cell, with one trap at the center of the 
unit cell, and the second one crossing the corner of the cell pe- 
riodically. (C) Illustration of Hund's first rule with repulsive 
contact interactions: the energy is lower as the particles be- 
come non-interacting for the same species, i.e. aligned spins 
(see text). 



tunneling can be tuned by varying the lattice depth Vq. 
With increasing Vq the atoms become more localized in 
the lattice sites, the band dispersion decreases and the 
shells in the individual traps are separated by increas- 
ingly large gaps. For sufficiently large lattice depth Vb, 
it is the shell structure of the the individual atom traps 
at the lattice sites that determines the physical behavior 
of the lattice. 

The mechanism leading to magnetic effects in the single 
quantum wells of the lattice is sketched in Figure 1(C). In 
the two-component contact-interacting fermion gas, due 
to the Pauli principle there is no mutual interaction be- 
tween atoms of the same species. Therefore, the trapped 
atoms in a degenerate shell can lower their interaction 
energy by maximizing the number of atoms of the same 
species - in other words, by aligning their spins. This 
mechanism in contact-interacting atomic systems leads 
to Hund's first rule 2 ^ and magnetism, in close similarity 
to long-range interacting electronic systems, as for exam- 
ple, quantum dots in semiconductor heterostructures 2 ^. 
Here, however, Hund's rule has a more dramatic effect 
since it completely removes the interaction between the 
same atom species. 

A Bose-Einstein condensate of a weakly interacting, di- 
lute gas of bosonic atoms is known to be well described 



by the Gross-Pitaevskii equation for the condensate 
wave function 2 *. Correspondingly, contact-interacting 
fermions may be approximated by a set of Kohn-Sham- 
like equations with a local effective potential. We note 
that in the dilute limit the exchange energy is treated 
exactly 25 . Periodic boundary conditions imply Bloch 
form for the orbitals, VVika-(r) = exp(ik-r)ii n i C(7 (r), where 
n labels the band, a = (J,, f) is the spin index and the 
wave vector k is confined into the first Brillouin zone. 
The periodic functions u„kcr(r) satisfy 

h 2 

- 7— (V + ik) 2 u„ kff (r) + [V opt (r) + gn a (r)]w nk(T (r) 
2m 

= enk CT Wnko-(r), a ^ a' (1) 

where m is the atom mass, n a is the density of atom 
species a and g is the interaction strength. For the lat- 
ter we have set g = 0.3 E^/k 2 in order to stay in the 
weak-interaction regime, where Eq. ^ is valid. This was 
confirmed by comparing the mean-field calculations for a 
single 2D harmonic trap with results obtained by exact 
diagonalization for N < 8 (see Y. Yu et ai, to be pub- 
lished) . For the Bloch wave vector we use a grid of 5 x 5 
points in the first Brillouin zone, and the key results were 
tested with up to 9 x 9 points. In the band-structure cal- 
culation, the functions M n ko-(r) are expanded in a basis of 
11 x 11 plane waves. The self-consistent iterations were 
started with antiferromagnetic and ferromagnetic initial 
potentials and small random perturbations were added 
to the initial guesses in order to avoid convergence into 
saddle points of the potential surface. 
The results of the calculations in the above scheme are 
displayed in Figure 2, which shows the systematics of the 
magnetism as a function of the shell filling of the single 
quantum wells in the lattice as a function of N and Vq- 
The calculated points and the different ground state con- 
figurations are indicated explicitly in the figure. Filled 
and open circles are the antiferromagnetic (AF) and fer- 
romagnetic (F) states, respectively. The non-magnetic 
region at smaller values of Vq in the phase diagram is 
shaded and the corresponding states are indicated by 
squares (NM). The different magnetic phases follow sys- 
tematically the filling of the single-trap shell structure. 
The bottom of the trapping potential can be approxi- 
mated by a harmonic potential with fru> = AVqEr. In 
two dimensions, this leads to the lowest closed shells at 
particle numbers A^ —2, 6, and 12 (which are frequently 
also called "magic numbers" in the literature, refcring 
back to the conventions in nuclear physics) 2 . At higher 
energies, however, the potential has a notch connecting 
the different lattice sites which imposes a square symme- 
try that breaks the three-fold degeneracy of the third, 
2sld, oscillator shell. Due to this notch, in our case the 
closed shells correspond to atom numbers A" —2, 6, 10 
and 12. In these cases, the spins in the single quantum 
wells are compensated, and the lattice is non-magnetic. 
However, in between these values, at N = 1, 4, 8 and 11, 
the shells are half-filled. Antiferromagnetic ordering of 
spins is observed as a gap is formed at the Fermi-level. 
For other atom numbers, the shells are only partially 
filled, and therefore the Fermi-levels reside in the mid- 
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FIG. 2: Magnetism in a square optical lattice, as a function 
of the lattice depth Vo (determining the tunneling and cou- 
pling between the single wells) in units of the recoil energy 
Er, and the particle number (determining the shell filling). 
The different non-magnetic (NM, blue squares), ferromagnetic 
(F, empty circles) and anti-ferromagnetic (AF, filled circles) 
phases are indicated. For closed shells ("magic numbers") 
at N = 2, 6 and 10, the lattice is non-magnetic. The anti- 
ferromagnetic ground state is seen at mid-shell only, i.e. at 
N = 1,4 and 8. In between, i.e. at the beginning and end of 
a shell, the lattice is ferromagnetic. The shaded area (green 
color) indicates the non-magnetic region at smaller values of 
Vb. 



die of a band. These configurations cannot open a gap. 
They favor ferromagnetism, as the exchange-like effect 
splits the spin levels. The lowest band is formed from 
the Is levels of the individual wells. It can be occupied 
by at most two atoms per lattice site. By assuming an 
antifcrromagnetic alignment of spins at half-filled band 
N = 1, a gap opens at the Fermi-level since the near- 
est similar neighbor is twice as far away as in the ferro- 
magnetic lattice where every lattice site is identical. For 
N = 3, 4, 5 and 6 the levels at the lp shell are occupied. 
At N = 4 the spin at the lp shell is maximized due to 
Hund's first rule. This mid-shell configuration favors an- 
tiferromagnetism where the energy is reduced by opening 
the Fermi-gap. At N = 3 the shell has only one atom 
in a p-shell and at N ~ 5 the shell is nearly filled. In 
both these cases the Fermi-energy lies on the band and 
ferromagnetism is found. 

At the third shell, 2sld, a shell closing at N = 10 with 
a non-magnetic phase, and mid-shell fillings for N = 8 
and N = 11 with antiferromagnetic phases are encoun- 
tered. This is attributed to the symmetry-breaking of 
the 2sld shell that resides close to the potential junction 
between different lattice sites. The d-state whose density 
lobes are not oriented towards the nearest neighbors is 
pushed higher in energy. As a result, the 2sld shell is 
split into two sub-shells, one consisting of 2s and Id lev- 
els and the other one formed by the higher-lying Id level. 
Consequently, we find mid-shell fillings at N = 8 and 11 
leading to antiferromagnetism. 

The splitting between the spin bands in the ferromag- 



FIG. 3: Integrated spin density over a single lattice site ("spin 
per site", Sg-^/h), as a function of increasing lattice depths 
Vq (in units of the recoil energy En). The different colors 
correspond to different particle numbers (shell fillings) of the 
single traps, as indicated in the inset. 



netic cases is of the order of 0.050 Er which correspond 
to approximately 0.048 ksTp in a gas of 40 K, where the 
Fermi-temperature is^li Tp — 330 nK . The typical en- 
ergy difference between non-magnetic and ferromagnetic 
states is 0.001 Er, corresponding to roughly 0.001 ksTp 
in the fermionic potassium gas. These energies corre- 
spond to temperatures which are about a factor of ten 
lower than the ones achieved experimentally. We men- 
tion here that an interaction-induced cooling mechanism 
was proposed^! in order to reach the antiferromagnetic 
phase at N = 1. 

The onset of magnetism as a function of Vq can be seen 
in Figure 3 which shows the " spin per site" , obtained by 
integrating the spin density over a single lattice site, 

%te = %\ W(v)-n\v)\dv. (2) 

* J site 

Generally, the magnetism sets on with increasing Vq as 
the number of atoms per lattice site becomes larger and, 
thus, the spatial extent of the highest occupied orbital in- 
creases. The transition occurs roughly at the same value 
of Vo when the Fermi-level resides at particular shells. 
For example, the transition occurs at Vq ~ 4.2 En for 
atom numbers N — 3, 4 and 5 as the lp levels are oc- 
cupied. For N = 11 atoms, antiferromagnetism sets on 
already around Vo = 5 Er. In this case, the highest 
occupied orbital (Id) has only a small overlap with the 
corresponding orbital at the neighboring site. Therefore, 
the optical lattice has to be shallow before the tunnel- 
ing between these levels is large for the magnetism to 
disappear. 

In summary, going beyond the single-band Hubbard 
model we found that two-component cold fermionic 
atoms in an optical lattice establish a rich magnetic phase 
diagram, and can show analogous effects to the electrons 
in magnetic solids. The magnetism in the optical lattice 
is determined by the number of atoms per lattice site, 
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and the depths of the lattice. If the lattice is shallow, for 
strong inter-site tunneling the lattice does not show any 
magnetism. For deeper lattices, where the tunneling is 
small, the total spin of the N atoms at the single sites 
is determined by the shell structure in the lattice min- 
ima. For closed shells (so-called "magic numbers") at 
N = 2, 6, 10 and 12, the single-wells are non-magnetic. 
For contact-interacting repulsive fermions, Hund's first 
rule applies in a particularly dramatic way, removing the 
interaction between the same atom species. Antiferro- 
magnetic ordering of the single-site spins was found when 
a gap at the Fermi level could open up. This was the case 
at mid-shell, with N = 1,4,8, or 11 particles in the sin- 
gle traps. Ferromagnetism occurred at the beginning and 
the end of a shell. 



FIG. 4: Spin density (n-f — njj of the ferromagnetic (F) and 
antiferromagnetic (AF) states for TV = 9 (Vo = 11-Ei?) and 
N — 8 (Vo — 13-Er), respectively. (The amplitude is renor- 
malized to maximum value at unity.) 
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